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Basic equations of nonequilibrium thermo field dynamics of dense quantum 
systems are presented. A formulation of nonequilibrium thermo field dynam- 
' ics has been performed using the nonequilibrium statistical operator method 

t^J- , by D.N.Zubarev. Hydrodynamic equations have been obtained in thermo field 

' representation. Two levels of the description of kinetics and hydrodynamics of 

a dense nuclear matter are considered. The first one is a quantum system with 
strongly coupled states, the second one is a quark-gluon plasma. Generalized 
q>^ , transfer equations of a consistent description of kinetics and hydrodynamics 

have been obtained, as well as limiting cases are considered. 
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1. Introduction 



The development of methods for the construction of kinetic and hydrodynamic 
equations in the theory of nonequilibrium processes for temperature quantum field 
systems is, in particular, important for the investigation of nonequilibrium properties 
of a quark-gluon plasma [ [l], @, 0, d, |] - one of the nuclear matter states which can 
be created at ultrarelativistic collisions of heavy nuclei [ |6|, [?|, |8|, |?J . The description of 
strongly nonequilibrium processes of a nuclear matter appears also in the laser ther- 



monuclear synthesis in systems D-D, D-T, D- 3 He, p-B [ |10], |TT], [12], |l3fl . In theoretical 
investigations of the laser thermonuclear synthesis there is a problem of description 
of a laser beam propagation and its absorption by the target (in particular, by the D- 
T mixture), electron and photon energy transport in an ionized magnetized plasma 
with the creation of a corona - an electron liquid (n e ~ 10 21 -j- 10 22 cm -3 ), and the 
creation of a core - a super dense ion liquid (n ~ 10 24 -=-10 26 cm -3 ). In target ablation 
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and implosion processes it is important to note transport coefficients, such as ther- 
mal conductivity, electron's conductivity, as well as dielectric properties, generation 
and evolution of spontaneous magnetic fields (with the inductivity of thousands 
of Tesla), nuclear synthesis mechanisms and energy transport by a-particles and 
neutrons which are products of synthesis reactions. In this connection, there is a 
series of important problems, when kinetic and hydrodynamic processes in a mag- 
netized super dense electron-ion plasma should be described consistently. Nonlinear 
hydrodynamic fluctuations of a magnetized degenerate electron liquid in a corona 
should be considered by taking into account both the kinetics of super heat electrons 
and nonlinear hydrodynamics of a super dense magnetized ion liquid, where nuclear 
synthesis reactions take place with the creation of high energy neutrons (with con- 
centration ~ 10 24 cm -3 ) and a-particles. The kinetics of these particles affects the 
synthesis processes and energy transport in the whole system. From the point of view 
of theoretical methods, in order to describe such strongly correlated nonequilibrium 
processes, one needs to construct the kinetics and hydrodynamics of a super dense 
high temperature electron-ion plasma and a nuclear matter at ultrarelativistic colli- 
sions of heavy nuclei. A nuclear matter - a stage of a quark-gluon plasma - appears 
after the compression of the target core of D-D, D-T plasma. As this takes place, its 
density increases by a factor of ten in the fourth degree and the distance between 
nuclons in the centre reaches ~ 10~ 13 cm. Such systems are examples of strong both 
long-range and short-range (nuclear) interactions. There is no small parameter for 
these systems (density, for example). Nonequilibrium processes have a strongly cor- 
related collective nature. That is why methods which are based on a one-particle 
description, in particular, on the basis of Boltzmann-like kinetic equations, cannot 
be used. In addition to high temperature dense quantum systems, there are Bose 
and Fermi systems at low temperatures with decisive many-particle dissipative cor- 
relations. Neither the linear response theory nor Boltzmann-like kinetic equations 
are sufficient for their description. 

Analysis of the problem of a description of kinetic processes in highly nonequi- 
librium and strongly coupled quantum systems on the basis of the nonequilibrium 
real-time Green functions technique [ [14], [1^, [L6| and the theory in terms of non- 



Markovian kinetic equations describing memory effects [ |17], |19|] was made in 



recent paper [ 2C] and then in monograph [21]. It is important to note that in [ 
20| the quantum kinetic equation for a dense and strongly coupled nonequilibrium 
system was obtained when the parameters of a shortened description included a one- 
particle Wigner distribution function and an average energy density. On the basis 
of this approach the quantum Enskog kinetic equation was obtained in [ ^TJ . This 
equation is the quantum analogue of the classical one within the revised Enskog the- 
ory [ ^3| . Problems of the construction of kinetic and hydrodynamic equations 
for highly nonequilibrium and strongly coupled quantum systems were considered 



based on the nonequilibrium thermo field dynamics in [ |24|, gg, |26, £7[ 13- ^ n 



particular, a generalized kinetic equation for the average value of the Klimontovich 



operator was obtained in [ |25[ with the help of the Kawasaki-Gunton projection 



operator method. The formalism of the nonequilibrium thermo field dynamics was 
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applied to the description of a hydrodynamic state of quantum field systems in 
paper [ ^6] . Generalized transport equations for nonequilibrium quantum systems, 
specifically for kinetic and hydrodynamic stages, were obtained in [ on the basis 
of the thermo field dynamics conception [ RTL using the nonequilibrium statis- 



tical operator method [ [H], |33|, |J. In this approach, similarly to [ |21| , the 
decisive role is that a set of the observed quantities is included in the description 
of the nonequilibrium process. For these quantities one finds generalized transport 
equations which should agree with nonequilibrium thermodynamics at controlling 
the local conservation laws for the particles-number density, momentum and energy. 
It gives substantial advantages over the nonequilibrium Green function technique [ 
T\ 0, which quite well describes excitation spectra, but practically does not 



describe nonequilibrium thermodynamics, and has problems with the local conser- 
vation laws control and the generalized transport coefficients calculation. 

In this paper we consider the kinetics and hydrodynamics of highly nonequi- 
librium and strongly coupled quantum systems using the nonequilibrium thermo 
field dynamics on the basis of the D.N.Zubarev nonequilibrium statistical operator 
method [ . Within this method we consider two different levels of a description of 
the kinetics and hydrodynamics of dense quantum nuclear systems: strongly coupled 
states and a quark-gluon plasma. The nonequilibrium thermo field dynamics on the 
basis of the nonequilibrium statistical operator method constitutes sections 2 to 7. A 
nonequilibrium thermo vacuum state vector is obtained here in view of equations for 
the generalized hydrodynamics of dense quantum systems. Transport equations of a 
consistent description of the kinetics and hydrodynamics in thermo field represen- 
tation are obtained in section 8. We mean that these equations are applied to dense 
quantum systems where strong coupled states can appear. This item implies, as one 
of the approaches, to investigate a nonequilibrium nuclear matter [ ^ ||. Another 
level of the description (rather a microscopic one) of nonequilibrium properties of 
dense quantum systems is considered in section 9. A consistent description of the 
kinetics and hydrodynamics is obtained here for a quark-gluon plasma. 

2. Thermo field dynamics formalism. Superoperators and state 
vectors in the Liouville thermo field space 

In this section we reformulate the nonequilibrium statistical mechanics of quan- 



tum systems using the thermo field dynamics formalism [ [31], [32 . 

Let us consider a quantum system of N interacting bosons or fermions. The 
Hamiltonian of this system is expressed via creation a z + and annihilation a t operators 
of the corresponding statistics: 

H = H(a + ,a). (2.1) 
Operators a z + , a l satisfy the commutation relations: 

[a t , a+] CT = Sij, [a t , a j ] a = [a+, a+] CT = 0, (2.2) 

where [A, B] a = AB — aBA, o = +1 for bosons and a = — 1 for fermions. 
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The nonequilibrium state of such a system is completely described by the non- 
equilibrium statistical operator g(t). This operator satisfies the quantum Liouville 
equation 

lg(t)-±[H,g(t)] = 0. (2.3) 

The nonequilibrium statistical operator g(t) allows us to calculate the average values 
of operators A 

(A) t = Sp(Ag(t)), (2.4) 

which can be observable quantities describing the nonequilibrium state of the system 
(for example, a hydrodynamic state is described by the average values of operators 
of particle number, momentum and energy densities). 

The main idea of thermo field dynamics [ [H], [32|] and its nonequilibrium formu- 
lation [ |35|, |36], [37| consists in doing the calculation of average values (|2.4j) with the 



help of the so-called nonequilibrium thermo vacuum state vector: 

(AY=((l\Ag(t))) = ((l\A\g(t)}), (2.5) 

where A is a superoperator which acts on state \g(t))). Nonequilibrium thermo vac- 
uum state vector \g(t))) satisfies the Schrodinger equation. Starting from equation 
r3~D, we obtain the relation 



3 W))) 



Of 

or, opening commutator, 

^\Q(t)))-^H\g(t))) = 0- (2-6) 

Here the "total" Hamiltonian H reads: 

H = H — H, (2.7) 

and it is known that ((1\H — 0; H — H(a + ,a), H = H*(a + ,a) are superoperators 
which consist of creation and annihilation superoperators without and with a tilde, 
and which represent the thermal Liouville space [ |39j . Superoperators H and H 
are defined in accordance with the relations: 

\Hg{t)))=H\g{t))), 

\g(t)H))=H\g(t))). { ' > 



Hence, it appears that at going from the quantum Liouville equation ( |2.3|) for 
nonequilibrium statistical operator g(t) to the Schrodinger equation ( |2.6| ) for non- 
equilibrium thermo vacuum state vector \g(t))), according to ( p.5|) , the number of 
creation and annihilation operators is doubled. Superoperators a+, , af, dj satisfy 
the same commutation relations as for operators a z + , ( |2.2| ): 

[d h a+] CT = [a h af] a = Sij, [d h dj) a = [a+, a+] CT = 0, (2.9) 
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[at, a,j\ a = [af, aj\ a = 0, [a h aA a = [af, aj\ a = 0. 
Annihilation superoperators a l , a l are defined in accordance with their action on 



the vacuum state - the supervacuum [ 38 



a,|00» =5,|00)) = 0, (2.10) 

where 1 00) ) = ||0)(0|)) is a supervacuum; and it is known that a t |0) = a l |0) = 0, 
and (0|fij = 0, i.e. a supervacuum 1 00) ) is an orthogonalized state of two vacuum 
states (0| and |0). Taking into account commutation relations ( |2.9| ), (|2.10|) , one can 



introduce unit vectors |1)) = IX/JOOD) anc ^ (UI = ((Xw 10(^11 m the following 
forms: 



|l»=exp £a+a+ 100)), 

I i J (2.11) 
«l| = «00|exp|Ea, a t }. 

With the help of these expressions one can find relations between the action of 
superoperators af, a - , af, dj 

a+|i» = ^ «i|a f = «i|a+«r. [ZAZ) 



In such a way, in the thermal field dynamics formalism [ [31], [32| the number of op- 
erators is doubled by introducing both without tilde and tildian operators A{a + , a), 
A(d + , a) for which the following properties take place: 

MA 2 = A X A 2 , A = A, 
c x A x + c 2 A 2 = c^Ai + c* 2 A 2 , / 2 13 n 

\A)) = A\1)), 
\A X A 2 )) = A l \A 2 )). 

Here * denotes a complex conjugation. Some detailed description of the properties 
of superoperators well thermal Liouville space is given in 

papers [0,0, 0,0. 

The nonequilibrium thermo vacuum state vector is normalized 

«l| ff (t)» = ((l|ff(t)|l» = l I (2.14) 

where g(t) is a nonequilibrium statistical superoperator. It depends on superopera- 
tors af, a L 

g(t) = g(a + ,d;t), (2.15) 

and, it is known that the corresponding tildian superoperator g(t) = g + (d + ,a;t) 
depends on superoperators af, 5^ . 
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3. Nonequilibrium statistical operator in thermo field space 

To solve the Schrodinger equation (|2.6Q a boundary condition should be given. 



Following the nonequilibrium statistical operator method [ [21], ^7], [34]], let us 
find a solution to this equation in a form, which depends on time via some set of 
observable quantities only. It means that this set is sufficient for the description of a 
nonequilibrium state of a system and does not depend on the initial moment of time. 
The solution to the Schrodinger equation, which satisfies the following boundary 
condition 



\g(t)))t=t = lf? q (*o)», (3-i; 



reads: 



\g(t))) = exp j(t - *o)^} |ft,(to)»- (3-2) 

We will consider times t 3> to, when the details of the initial state become inessential. 
To avoid the dependence on t , let us average the solution (|3.2|) on the initial time 
moment in the range between i and t and make the limiting transition t —t — > — oo. 



We will obtain [|7|: 



o 



\g(t))) = e J dt' e £t e--^\ Qq (t + t'))), (3.3) 

— oo 

where e —>■ +0 after the thermodynamic limiting transition. Solution ( |3.3| ), as it can 
be shown by its direct differentiation with respect to time t, satisfies the Schrodinger 
equation with a small source in the right-hand side: 

I - k R ) m)) = - £ (^ (t))) ~ (3 - 4) 

This source selects retarded solutions which correspond to a shortened description of 
the nonequilibrium state of a system, \e q (t))) is a thermo vacuum quasiequilibrium 
state vector 

Ifiq(*)» = fiq(*)|l»- (3-5) 

Similarly to ( [2.14| ), it is normalized by the rule 

((l|^ q (t))) = ((l|^(t)|l)) = l, (3.6) 

where g q (t) is a quasiequilibrium statistical superoperator. The quasiequilibrium 
thermo vacuum state vector of a system is introduced in the following way. Let 
(PnY = ((l|Pn| £?(£))) be a set of observable quantities which describe the nonequilib- 
rium state of a system. p n are operators which consist of the creation and annihila- 
tion operators defined in fl2.2|) . Quasiequilibrium statistical operator g q (t) is defined 
from the condition of informational entropy S^f extremum (maximum) at additional 
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conditions of prescribing the average values (p n Y and conservation of normalization 
condition Q [0,0: 



g q (t) = exp - j 

d>(t) = InSpexp j-^F*(t)p n J, 



(3.7) 



where $(t) is the Massieu-Planck functional. A summation on n can designate a sum 
with respect to the wave- vector k, the kind of particles and a line of other quantum 
numbers, a spin for example. Parameters F n (t) are defined from the conditions of 
self-consistency: 

(p n Y = (p n )% (...>i = Sp(... fiq (*)). (3.8) 

According to ( |2.5|) , let us write these conditions of self-consistency in the following 
form: 

= «l|ft»|flq(f)»- (3-9) 

Taking into account behaviours ( [2. 131 ), we have: 

lfiq(*)» = fti(*)|i» = ef(*)|i». (3-io) 



where 



QS) =exp|-^(*)-Z)^:(*)p»J 

g+(t) = exp j-$(t) - J2 F n{t)Pn\ 



(3.11) 



are quasiequilibrium statistical superoperators which contain superoperators p n and 
p n , correspondingly: 

Pn = P n (a + ,S). 

If self-consistency condition ( |3.9|) realizes, we shall have the following relations (at 
fixed corresponding parameters): 

HHf) -((l\p n \ 6(l (t))) = ((MPn\Q(t))), 



SFjt) 



(3.13) 

" <1>(/) -((l\Pn\Q q (t))) = ((MPn\e(t)))- 
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Relations Q3.13|) show that parameters F*(t), F n (t) are conjugated to averages 



((l\p n \g(t))) and ((l\p n \g(t))), correspondingly. On the other hand, with the help of 
\g q (t))) and self-consistency conditions ( |3.9|) we can define the entropy of the system 
state: 

S(t) = g q (t))g q (t))) = $(t) + £ i^(t)«l|pn|e(*)». (3.14) 

n 

The physical meaning of parameters F*(t) can be obtained now on the basis of the 
previous relation: 

K(t) = imh- (3 - 15) 

4. Projection operators in thermo field space 

Now the auxiliary quasiequilibrium thermo vacuum state vector \g q (t))) is de- 
fined. Let us represent solution ( |3.3|) of the Schrodinger equation ( |2.6| ) in a form 
which is more convenient for the construction of transport equations for averages 
((l\p n \g(t))) . We shall start from the Schrodinger equation with a small source (|3.4| ). 
Let us rebuild this equation by introducing A\g(t))) = \g(t))) — \g q (t))): 



The calculation of time derivation of \g q (t))) in the right-hand side of equation ( |4.1|) 
is equivalent to the introduction of the Kawasaki- Gunton projection operator p q (t) 
[ |27j in thermo field representation: 

= P q {t)^H\g{t))), (4-2) 
p q (t)(\...))) = \g q (t)}}+ (4-3) 

Projection operator p q (t) acts on state vectors | . . .)) only and has all the operator 
properties: 

P*(t)W))) = l<?q(*)», 

pMe^t'))) = |f? q (*))>, 

Pq(*)Pq(0 = P$)- 

Taking into account condition p q (t)rgH A\g(t))) = 0, one may rewrite equation 
Jj), after simple reductions, in a form: 



(| " (l " P*®)fiH + £ ) = (l - p q (t)) ^H\g q (t)}}. 



(4-4) 
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The formal solution to this equation reads: 



A|f?(*)» 



t 

J dt' e^^T(t,t')[l-p q (t'))^H\g q (t'))), 



or 



(*)» = Mt)))+ j dt'e^'- t )T(t ) t')(l-p q (t , ))^b q (t , )» ) (4-5) 



where 



T(t,t') = ex P+ y*dt' (l - p q (f))^#} 



(4.6) 



is an evolution operator with projection consideration, and exp + is an ordered ex- 
ponent. Then, let us consider expression ^1 — p q (t')j ^H\g q (t))) in the right-hand 

side of ( |4.1|) . The action of j^H and ^1 — p q {f)\ on \g q (t))) can be represented in 



the form 

f 



-Pq(0)^lffq(*)» = E J? »( t ) /^^(l-pCO)^^))), (4-7) 



where p n and p(t) read: 

Pn = -Tz[H,p n ], 
in 

P(t)p = (WMt))) + E ^|^gft))) (^ - ((MPnlsit)}) 



(4.8) 
(4.9) 



Here p(t) is a generalized Mori projection operator in thermo field representation. 
It acts on operators and has the following properties: 



(4.10) 



P(t)p n = Pn, 
p(t)p(t>) = p(t). 

Let us substitute now (|4.7D into (|4.5| ) and, as a result, we will obtain an expression 
for the nonequilibrium thermo vacuum state of a system: 



t 

!<?(*)» = lflq(*)» + £ / dt'e^T(M') 



fdrgl(l/)J n (l/)g^(l/)))F:(l/) 



Here 



J n (t) = (l - p(t))p n 
are generalized flows. 



(4.11) 
(4.12) 
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Let us obtain now transport equations for averages ((l\p n \g(t))) in thermo field 
representation with the help of nonequilibrium thermo vacuum state vector \g{t))) 
( |4.11|) . To achieve this we will use the equality 

|«l|Pn|^)» = = ((MkMt))) + ((Jn(t)\g(t))}. (4.13) 

By making use of \g{t))) in Q4.11| ) in averaging the last term, we obtain transport 
equations for ((l\p n \g(t))): 

d 



Of 



((llAil *?(*)» = (WMt))) 



(4.14) 



£ / dt ' e*'-* ((j n (t)T(t,t') 



fdr gl(tVnit'H~ T (t')))F:,(t 



where p n = — ^Hp n . Relations (|4.14|) are treated as a general form of transport 
equations for average values of a shortened description. These equations can be 
applied to completely actual problems. 

In the case of weakly nonequilibrium processes, the generalized transport equa- 
tions Q4.14D are reduced appreciably. We shall consider this case in the next section. 



5. Transport equations in linear approximation 



Let us consider the nonequilibrium state of a quantum field system near equilib- 
rium. In this connection let us suppose that average values {{l\p n \g(t))) of variables 
for a shortened description and their conjugated parameters F*(t) differ in their 
equilibrium meanings slightly. In such a case, one can expand the quasiequilibrium 
thermo vacuum state vector \g q (to))) ( |3.5| ), ( |3.11| ) into a series on deviations of 
parameters F*(t) from their equilibrium values F n (0) and restrict a linear approxi- 
mation only: 



\g q (t))) = \Qo(F n (0))))-J2 6F n^ 



fdr ^n(0)K^(F n (0)) 



(5.1) 



where \go{F n (0)))) is an equilibrium thermo vacuum state vector which depends on 
equilibrium values of F n (0) parameters (local inversed temperature f3 and chemical 
potential /x), 5F*(t) = F*(t) — F*(0). With the help of self-consistency conditions 
( P-9|) and taking into account Q5.1]) , let us define parameters SF*(t): 

(5.2) 



SF:(t)^-J2((p\p))nL((MSPm\e(t))), 



where 5p m = p m — ((l\p m \go(F n (0)))), {{p\p))^ n are elements of the inverse matrix 
of ((p\p)). Elements of the matrix {(p\p)) are equilibrium correlation functions in 
thermo field representation 



{(Pn\Pm))o = ({Pn 



T QoPmOo T 



(5.3) 
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Substitution of (|5.3|) into ( |5.1|) for \g q (t))) results in 



|0q 



(t))) = \qo(F(0)))) +Y,(( 1 \5Pm\Q(t)))((p\p))^ n 



fdr ^(F(O)K^^(F(O)) 



(5.4) 



To calculate nonequilibrium thermo vacuum state vector \g(t))) (|4.11 ) in linear 
approximation ( |5.4|) , let us rewrite transport equations ( |4.14|) in another form: 

^«l|*Pn|<?(*)» = J2n nm ((l\6p n \g(t))) fdt> e £ (*'-')^ m (t,t')((l|^ n .|^ / )}), 



where 



fdr elmel" }) ((*)),;! 



are equilibrium quantum correlation functions, and 



ip 7 



l (t,t') = E(( J " 



T Q (t,t') 



fdr elJtifr 7 )) ((p\p))^ 



(5.5) 



(5.6) 



(5.7) 



are transport cores in thermo field representation, which describe dissipative pro- 
cesses in the weakly nonequilibrium state of a quantum system, J n is defined simi- 
larly to flOp : 

J n = (1 - p )pn; 



T (t) = expU(l-po)-H 

is the time evolution operator with taking into account the projecting where po is 
the Mori projection operator with the following structure: 



Po A=((l\A\g )) + J2mPnMp\p))- 1 m 



rnPrn- 



(5i 



Projection operator po satisfies conditions (|4.10[) . As it can be shown [ P7H , for trans- 
port equations (|5.5| ) correspond to equations for quantum time correlation functions 
in thermo field representation: 



0_ 

Of 



® nm {t) = ((Pn(t)\p m (0)g )), (5.9) 

t 

®nm(t) = J2 n nl®lm(t)-J2 [ dt ' ^"'W^ t')®lm(t'), (5.10) 

7 7 J 



where the time evolution of superoperators p n (t) in the Heisenberg representation 
reads: 



1 .-. 



p n (t) = exp <^ -TrHt ) p n exp <^ —Ht 



ih 



l 



ih 



(5.11) 
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As it is indicated by ( |5.5|) , (|5.10| ), transport equations for values ((l\p n \g(t))) and 



the corresponding equations for time correlation functions $ nm (t) in linear approx- 
imation (|5.4| ) are closed. 

If spectral magnitudes $ nm (u;) for quantum time correlation functions $„ m (t) 
are defined as follows: 



$nm(w) 

Ztt 



oo 

l - J dte iu>t $ nm (t), (5.12) 



then, according to the definition of Green functions [ |33[ and the condition that 
Qnm{^) is a real function (u; is a real number), one finds imaginary parts of spectral 
magnitudes of the corresponding retarded, advanced and causal Green functions 

5mCM = (e^ - a) $„ m (^), (5.13) 

^mffjw) = ~\ (e^ + cr) $ nm (a;), (5.14) 

where functions G T nm (t,t'), G^ m {t,t'), G^ m (t,t') in time representation constitute, 
correspondingly, retarded, advanced and causal Green functions in thermo field rep- 
resentation: 

G^(M') = ±^0(t-t')((l\[ Pn (t),p m (t')]M) = (5.15) 

±^e(t-t')((i\[p n (t),p m (t%\go)) = 

±Tj.e(t - t')((l\p n (t)p m (t') - ap m (t')p n (t)\g )), 

G c nm (t,t') = ^((l\Tp n {t)p m (t')g )) (5.16) 

= h(t - t'){{l\p n {t)p m {t')\ Qo )) + - t')((l|^(*)Pm(OkQ)), 
m in 

where 6{t) = | q' ^ < q' * s a un ^ s ^ e P function. As it is seen from (|5.15|) and ( |5.16|) , 

retarded (advanced) and causal Green functions are defined via time correlation 
functions ((l\p n (t)p m (t')\go)) , {{l\pn(t)p m (t')\go)), which satisfy conditions ( |5.10|) . 



For completeness, let us write dispersion relations which connect the imaginary parts 
SmG^(w), QmG^ m (w) (|5U3|) , Q514D and the real parts 9foG£* (w), UzG c nm {uj) 



of spectral functions of the corresponding Green functions [ [H], |33| : 



oo 

SeG^(w) = ±-p [ du/ ^ mGT '^'\ (5.17) 

71 J 10' — U 

— oo 
oo 

^tG c nm (oo) = -p / dtu — -^—^r, • 

nmK ' tt oj' -uj +cr 
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It should be pointed out that Green functions in thermo field representation, which 
are calculated with the help of equilibrium thermo vacuum state vector \qo(/3))), were 



investigated in many papers [ 32, B8L HQ. 41 1. Specifically, the diagram technique for 



their calculation, which generalizes the Feynman method, was developed in [ f|!J . 

An application of the general structure of nonequilibrium thermo field dynamics 
on the basis of the nonequilibrium statistical operator method [ |2l], [27| will be 
considered in sections 7 and 8. We will obtain equations of hydrodynamics and 
transport equations of a consistent description of the kinetics and hydrodynamics 
for dense quantum systems in thermo field representation. While investigating such 
systems, one of important problems is the calculation of transport cores (or transport 
coefficients) both for weakly and strongly nonequilibrium systems. 

The problem is that quasiequilibrium thermo vacuum state vector \g q (t))), ( |3.10|) 



or fl5.4p, in each case is not a vacuum state for superoperators d, d + , d, d + . The 
subject of this question consists in the construction of dynamical reflection of su- 
peroperators d, d + , d, d + into superoperators for "quasiparticles" for which the 
quasiequilibrium thermo vacuum state is a vacuum state. One method of construct- 
ing such superoperators is considered in the next section. 

6. Creation and annihilation superoperators of "quasiparticles" 
for the quasiequilibrium thermo vacuum state 



In our further consideration we assume that Hamiltonian H (|2.1|) of a system 
can be represented in the form: 

H = H + H inU (6.1) 

where H int contains a small parameter. This small parameter can be used for the 
construction of perturbation theory series. At the same time Ho is a non-perturbed 
part of the Hamiltonian (|6.1| ). It depends on creation and annihilation operators a z + , 
a l bilinearly. According to the nonequilibrium thermo field dynamics formalism, 
Hamiltonian H of a system reads: 

H = H +H int , 

H = Ho -Ho, (6.2) 
H m t = H mt —H mt . 

Hq depends on superoperators d + , d, d + , d bilinearly. In such a case, to construct 
the perturbation theory for operators it is convenient to use the Heisenberg repre- 
sentation on the non-perturbed part of Hamiltonian Hq: 

A(t) = exp \~H t \Aexp I ^H t\ . (6.3) 



ih J [ ih 

It is known that A(t) satisfies the Heisenberg equation 

j t A(t) = -i [H int (t),A(t)] , H int {t) = exp |-^#o*j ^exp l^H t\. (6.4) 



13 



.V.Tokarchuk, T.Arimitsu, A.E.Kobryn 



Superoperators a + , a, a + , a in the Heisenberg representation read: 



a(t)= e-i^aei^ *, a+(t) = e-^ Bot a + eTn Bot , 
5(f) = e-TK not aeTh Rot , a+(t)= e-Th Sot d + e^ n<}t . 



(6.5) 



They satisfy the "classical" commutation relations: 

[a(t), a + (t)] CT = [5(t), a + (t)] CT = 1, [a(t), a(t)] ff = [a + (t), a + {t)] a = 0, (6.6) 

[a(0, a(t)] a = [a + (t), a + (t)] CT = 0, [a(t), a{t)] a = [a + (t), a + (t)} a = 0. 

Let us assume that the quasiequilibrium thermo vacuum state vector describes 
the initial state of a quantum system by the non-perturbed part of Hamiltonian Hq. 
\g q (to))) is not a vacuum state for annihilation superoperators a(t), a(t), i.e. 

fi(*)lfti(*o)» = f{t-t )a + {t)\g q {t ))), 

a(t)\g q (t ))) = (Tf(t-t )a + (t)\g q (t ))), (6.7) 
<<l|a + (t) = a((l\~a+(t), 

where function f(t — to) will be defined below. Nevertheless, the linear combination 
of superoperators a + , a, a + , a allows us to define new creation and annihilation 
superoperators 7 + , and 7 + (t), j(t) [ [3^, |36|, |3j 



7(t) = Q 1/2 (t - t )[a(t) - f(t - t )ti + (t)], (6.8) 
7 + (t) = g 1/2 (t - t )[h + (t) - aa(t)\. (6.9) 



So, taking into account ( |6.7|) , the action of new operators on state vectors reads: 
7(*)l*?q(*o)» = 0, ((l| 7 + (t) = 0, 

7(*)|fi q («o)» = 0, ((1|7 + W = 0. lb ' iUj 

Superoperators 7 + , j(t) and 7 + (t), 7(t) satisfy the "canonical" commutation rela- 
tions: 

[7(t),7 + (t)] CT = 1, 

[7(t),7 + W]. = 1, [7(*),7(*)]a = [7 + (t),7 + W]. = 0. (6.11) 

A connecting expression between multiplier Q(t — to) an d function f(t — t ) may be 
found on the basis of relations ( p.8|) , ( |6.9|) and ( |6.11|) : 

Q(t-to) = [l-(T/(t-t )] _1 . (6.12) 



But, to define f(t — t ) function let us use the second equality in ( p.7|) . It aids to 
obtain the following: 

({l\a{t)a{t)\g q {t ))) = af{t - t )((l\a(t)a + (t)\g q (t ))). (6.13) 

And, using the third equality in fl6.7|) and ( |6.13| ) we arrive at 

((l\~a + (t)~a(t)\g q (t )}} = f(t - t )((l\a(t)a + (t)\g q (t )}}. (6.14) 
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As far as 

n{t - t ) = ((l\a + (t)a(t)\g q (t ))) = ((l\a + (t)a(t)\g q (t ))) (6.15) 
is the average particle number, then, using in the right-hand side of (|6.14j) commu- 



tation relation a(t)a + (t) — aa + (t)a(t) = 1, we obtain a linkage between n(t — to) and 
function f(t — t ), and vice versa: 

n(t-t ) = f(t-to)(l + <m(t-to)), 

f(t - t ) = n(t - t )/(l + an(t - t )). (6.16) 

Finally, if to substitute fl6.16j ) into (|6.12|) , one finds a linkage of normalized multiplier 
Q(t — to) and average particle number n{t — to): 

Q(t-t ) = l + an(t-t ). (6.17) 

Now, taking into account fl6.16|) and (|6.17|) , relations ( |6.8D , ( lof) for superoperators 
7 + , 7(t), and 7 + (t), j(t) read: 

n(t - t ) 



l(t) 
l + (t) 



1 + an(t — to 
1 + anit — to] 



a(t) — , 

_ w 1 + an(t - t 
a + (t) - aa(t)}. 



-a + {t) 



(6.18) 



Inversed transformations for superoperators a + , a(t) and a + (t), a(t) can be easily 
obtained from (|6.18| ): 

n(t - 1 ) 



a{t) = 1 + onit - to) 



lit) 



1 + unit — to) 
a+(t) = [l+an(t-t )V tf + (t) + a^{t)}. 



Tit) 



(6.19) 



j(t), 7 + (t), lit), 7 + (t) can be defined as annihilation and creation superoperators 
of quasiparticles for which quasiequilibrium thermo vacuum state \g q (t))) is a vac- 
uum state. These superoperators are functions of thermodynamic parameters F n (t) 
which describe the quasiequilibrium state of a system. Relations ( |6.18|) , (|6.19|) are 
a dynamical reflection between a(t), a + (t), a(t), a + (t) and *f(t), 7 + (t), lit), 7 + (t). 
While calculating transport cores in transport equations ( |4.14| ) and ( |5.7|) , or con- 
nected with these quantities Green functions, superoperators a(t), a + (t), a(t), a + {t) 
are multiplied according to the Wick theorem about normal products :a + (t)a(t):, 
\a + {t)a{t)\. Coming now to superoperators j(t), 7 + (£), l(t), l + (t) with the help 
of inversed transformations ( |6.19| ), one can obtain for the calculation of transport 
cores or the corresponding Green functions, a generalization of the Wick formulae 
and Feynman diagram technique on a nonequilibrium case. It is one of important 
features of the nonequilibrium thermo field dynamics. It consists in the possibility 
of constructing the quasiequilibrium thermo vacuum state vector as a vacuum state 
vector for the solution to the Schrodinger equation ( |2.6|) of a quantum field system. 

In the following section we consider equations of thermo field hydrodynamics 
of quantum systems of both strongly and weakly nonequilibrium states using the 



base transport equations Ql.l4| ). We obtain expressions for generalized transport 
coefficients of viscosity and thermal conductivity in thermo field representation. 
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7. Generalized hydrodynamic equations in thermo field represen- 
tation 



In the case of description of the nonequilibrium state of a quantum system the 
following quantities can be chosen as parameters of a shortened description: they 
are density operators of particle number n k , momentum p k and energy E k . For such 
a set of values quasiequilibrium statistical operator g q (t) reads [[21], |33|, |34| : 



g q (t) = exp \ -<fr(i) 



P- h (t)E(t) - (f3v)k{t)P k -P[H- jV 2 )\ k (t)h k 



m. 



(7.1) 



where parameters //_*.(£), V- k (t), /3-k{t) are defined from the self-consistency con- 
ditions 



((l|n fc |^(t))) =((l|n fc |^))), 

«l|p fc |eq(*)» = «l|Pfc 

{{l\E k \g q (t)))={{l\E k \g(t))), 



(7.2) 



and mean a chemical potential, average hydrodynamical velocity and a local value 
of inversed temperature, respectively [ |ST|. Expression ( |7.1| ) aids to obtain a set of 
equations of generalized hydrodynamics in thermo field representation. Taking into 
account ( f4.14p , we obtain: 



^«l|n*|ff(f)» = ((l\fi k \g q (t))), 
^((l\p k \g(t))) = ((l\ Pk \g q (t))) ~ 



E / dt' (/j p (k;t) 

° — oo 

t 

E / dt'e^'-)^J p (fc;t) 



r(t,f) 



T(M') 



Jdr ^(0 Jp(K *'K~ T (*') )> Wv) g V) + 



^((l\E k \g(t))) = ((l\E k \g q (t))) 



E / dt' (/j E (k;t) 



-oo 

f 



E / di ' e e <*V/j*(M) 



T(M') 



T(f,f) 



(7.3) 
(7.4) 



Jdr gl(t')J E (g;t')g q -^t') » /3 g (i'), 



(7.5) 



Jdr^(t')J P (g;t')^ r (t') » (Pv) g (t') + 



Jdr e ;(OJ B (g;Ofii" T (0)> ^ 
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where 



hk = n k (a£,o fc ) , 

Pk =Pk(at,a k ), J P (k;t) = (l — p(t)jp k , 
E k = E k {a+,a k ), J E (k;t) = (l-p(t))E k , 



are generalized flows. But the Mori projection operator p(t) at the description of 
the hydrodynamic state reads: 



P (f)(...) = «i|... Yl 

a={n,p,E} k 



g((l|---g q (t) 



a h -((l\a k \g(t)))) 



(7.6) 



Transport cores 

J p {k;t) T(t,t' 

J p {k-t) T(t,t' 
J E {k; t) 

J E (k;t) 



T(t,t>) 
T(t,t') 



Jdr Ql(t')J p (g;t>)^(t') 
o 

Jdr Q\{t')Mg;,t')Q\- T {t') 
o 

Jdr ^(0^(KOei _T (0 

o 



(7.7) 



in (|7.4|), ( |7.5[ ) are calculated with the help of quasiequilibrium thermo vacuum state 
vector ( [7.1|) . 

For the description of nonequilibrium hydrodynamic state of a quantum field 
system near equilibrium, the set of equations of generalized hydrodynamics ( |7.3|) - 
( |7.5|) according to ( ^Il) -( ^5|) becomes closed: 







-((l\5n k \ Q (t)}}=J2^nAk,g)((l\SpM t )))^ 
g 

g 

J2^ph(k,g)((l\Sh g \g(t)))- 
g 

t 

J dt' e^'^k : V (k,g;t,t') : g((l\5p g \g(t'))) 



(7.8) 
(7.9) 



-oo 

t 



1 dt ' e£(t '~ t)k ■■ U(*.k*.0 • g((i\shMt'))), 
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d 



(7.10) 



t 

E / dt ' e£[t '~ t)k ■ w*.«*.o : *<wje(0» - 



-oo 

/ 



£ / dt ' e^k - X(k,g,t,t') ■ g((l\5h g \g(t')))- 



In these relations all quantities Q np (k,g), Q pn (k,g), Q p h(k,g) and Qh P (k,g) are 
defined by the relation (|5.6| ) and are equilibrium quantum correlation functions; 
rj(k,g;t,t'), X(k,g;t,t'), £ ph {k, g;t,t'), £ hp {k, g;t,t') are, correspondingly, time-de- 
pendent generalized transport coefficients of viscosity, thermal conductivity and 
cross transport coefficients. They are defined as follows: 



k : rj(k,g;t,t') : g 
k : £ ph (k,g;t,t') ■ k 
k ■ £h P (k, g t, t') : g 

k-X(k,gt,t')-g 



J p (k)T (t,t') 
J p (k)T (t,t') 
J h (k)T (t,t') 
J h (k)T {t,t') 



Jdr QlJ P {g) Q l- T ))X pp (g 



fdr QlJ h {g)g 



l-T 



fdr glJ p {g)g 



l-T 



X pp {,§\ 



Jdr g^Mg^^X^g) 



(7.11) 
(7.12) 
(7.13) 
(7.14) 



where 



Jp(fe) = (1 - Po)P k , 
J h {k) = (1 - p )h k , 



h k = E h - XEn(k)X n n(k)n k , 



(7.15) 



and X PP (k), Xhh(k), X £«.(&), X nn (fe) are equilibrium quantum correlation functions: 



Xppig) 

Xnn {.&) 



fdr g T p_ g g 



l-T 



fdr g T n- s gl 



l-T 



Xhh{S) 

X En {g) 



Jdr g T Q h- g g 



l-T 



fdr g T n_ g gl 



l-T 



In ( |7.15| ) the Mori projection operator p is defined in accordance with ( |5.8| ) and is 
built on operators n k , p k , h k in thermo field representation. Operator h k appears 
due to the inclusion into \g q (t))) ( ff.l|) parameters f3 k {t) and fi k (t) with the help of 
self-consistency conditions (|7.2|) in the linear approximation (see ( |5.1| )- (|5~4j )) and 
has the meaning of generalized enthalpy. On the basis of hydrodynamic equations 
( [7.8| )-( |7TT0"D one can investigate quantum correlation functions like "density-density" 
$ nn (fc;t), "flow-flow" $p P (fc;t), "enthalpy-enthalpy" &hh(k;t) and the correspond- 
ing Green functions for specific quantum field systems. 
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A separate and very important problem is the calculation of transport coefficients 
( [7. 1 1|) — ( [7. 14j ) . Consideration of nonequilibrium thermodynamics in a field theory on 
the basis of nonequilibrium statistical operator [ [33] and approximate calculation 
of viscosity and thermal conductivity coefficients for the (£> 4 field model were done 
in paper [ |43fl. These calculations were carried out using the Green-Kubo formulae 
which connect transport coefficients of viscosity and thermal conductivity with the 
corresponding Green functions. The last ones are built on the stress tensor and 
energy flow operators. In their turn, the Green functions were calculated using the 



Dzyaloshinski diagram method [ 03 



8. Transport equations of a consistent description of the 
kinetics and hydrodynamics of dense quantum systems 

In the studies of nonequilibrium states of quantum field systems, such as a nu- 
clear matter [ ||, |9|, there arises a problem of taking into consideration coupled 
states. Kinetic and hydrodynamic processes in a hot, compressed nuclear matter, 
which appears after ultrarelativistic collisions of heavy nuclei or laser thermonuclear 
synthesis, are mutually connected and we should consider coupled states between 
nuclons. This is of great importance for the analysis and correlation of final reac- 
tion products. Obviously, a nuclon interaction investigation based on a quark-gluon 
plasma is a sequential microscopic approach to the dynamical description of reac- 
tions in a nuclear matter. The problem of a consistent description of the kinetics 
and hydrodynamics of a dense quark-gluon plasma is considered in section 9. For 
the description of kinetic processes in a nuclear matter on the level of model inter- 
actions, the Vlasov-Ueling-Uhlenbeck kinetic equation is used. This equation is used 
mainly in the case of low densities. The problems of a dense quark-gluon matter 



were discussed in detail in [ §, [| g(| [45], [4g, |47fl. 

We will consider a quantum field system in which coupled states can appear 
between the particles. Let us introduce annihilation and creation operators of a 
coupled state (Aa) with A-particle: 



a Aa 



[p) = Yl ^a P (l, • • • > A)a(l) . . . a(A), 



<M = £^ ap (l,...,A)a + (l)...a + (A), 



1 A 



where ^ Qp (l, . . . , A) is a self-function of the A-particle coupled state, a denotes 
internal quantum numbers (spin, etc.), p is a particle momentum, the sum covers 
the particles. Annihilation and creation operators a(j) and a + {j) satisfy the following 
commutation relations: 

[a(l), a + (j)] a = S hj , [a(0, a(j)] a = [a + (l), a + (j)] a = 0, (8.2) 

where a-commutator is determined by [a, b] a = ab — aba with a = ±1: +1 for bosons 
and —1 for fermions. 
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The Hamiltonian of such a system can be written in the form: 

^£/S£<( P -fk„( P+ ?)+ (8.3) 



(2irh) 6 2m A Aa V 2 J Aa V 2. 
Iv^v^ f dpdp'dq + ( q-p'\ A , . / q~p'\ 



where Vab(q) is interaction energy between A- and -B-particle coupled states, q is 
a wavevector. Annihilation and creation operators a Aa (p) and a Aa (p) satisfy the 
following commutation relations: 



Aa(P)i a B/3(P% = 5A,B5 a ,p5(p - p'), 
Aa(P), a B?(P% = [ a tM: a Bp(P% = °" 



(8.4) 



hspiq) in ( p.3|) is a Fourier transform of the -B-particle density operator: 
~ f \ f dp + 

HbM = J ppVfVf 

As parameters of a shortened description for the consistent description of the 
kinetics and hydrodynamics of a system, where coupled states between the particles 
can appear, let us choose nonequilibrium distribution functions of A-particle coupled 
states in thermo field representation 

((Mn J ^(r,p)\g(t))) = f^(r,p]t) = f Aa {x- i t), x = {r,p}, (8.5) 

here /^(a?; t) is a Wigner function of the A-particle coupled state where 

Tg^g ^^'^la (p - |) «4a (p + |) (8-6) 

is the Klimontovich density operator; and the average value of the total energy 
density operator 

((l\H(r)\g(t)}) = ((l\H(r)g(t))). (8.7) 

By this J dr H(r) = H, H{r) is a superoperator of the total energy density which 
is constructed on annihilation and creation superoperators a Aa (p) and a Aa (p). The 
latter satisfy commutation relations (|8.4j ). Following [27], one can rewrite quasiequi- 
librium statistical operator g q (t), \g q (t))) = g q (t)\l)) for the mentioned parameters 
of a shortened description in the form: 



g q (t) = exp < 



-&(t)-J dr/3(r;i) (#(r)-^ j H Aa (x; t)h Aa {x) 
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where Lagrange multipliers /3(r;t) and fi Aa (x;t) can be found from the self-consis- 
tency conditions, correspondingly: 



«l|ff(r)|e(t)» = ((l\H(r)\g q (t))), 
((l\n Aa (x)\g(t))) = ((l\n Aa (x)\g q (t))), 



(8.9) 
(8.10) 



is the Massieu-Planck functional and it can be defined from the normalization 
condition (|3.6|): 



= ln(( 1 



exp | -f dr P(r; t) \H(r) - £ J 



dp 



(2tvH) 



(8.11) 



Using now the general structure of nonequilibrium thermo field dynamics ( |4.1| )- 
( (4.14|) , one can obtain a set of generalized transport equations for A-particle Wigner 
distribution functions and the average interaction energy: 

d 



dt 



((Mn Aa (x)\g(t))) = ((l\h Aa (x)\g q (t))) + 



J dr' J dt' tfV-Q<ffe(x, r'- t, t')(3{r'- t') + 

— oo 

t 



1.12) 



0_ 

Of 



({l\H(r)\g(t))) = ((l\H(r)\g q (t))) 



J dr' J dt' e^'-^HHir, r'; t, t')(3(r'; t') + 



B,0 J 



where x' = {r',p'}, dx' = (2irh) 3 dr' dp'. Here 



<p^(x,r';t,t' 



<f H p n {r ,x;t,t 



i 



J n Jx,t)T(t,t') Jdr g T q (t')J H (r';t')g q -^t') )) , (f 



i 



"B0 



<p HH (r,r';t,t) = (( 1 



J H (r,t)T(t,t') Jdr g T q (t')J H (r';t')g q - T (t') 



.13) 



^(x,x';t,t') = {{lJ n Jx,t)T(t,t') Jdr gl(t')J n Jx';t')gl- T (t'))),(8.U) 

115) 

J H (r,t)T(t,t') Jdr gl(t')J n Jx';t')g q - T (t') )) , (8.16) 

5.17) 
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are generalized transport cores which describe dissipative processes. In these formu- 
lae 



Jnir-t) 

J n . Aa (r,p;t) = 

are generalized flows, 
H(r) - 



l-p(t'))H(r), 
1 - P(t'))h Aa (x) 



■—[H, n Aa (x)), 



p(t) is a generalized Mori projection operator in thermo field representation. It acts 
on operators 



P (t)P = ((\P\qM) + 
8((l\PMt))) 



119) 



dr 



E 

A, a 



5((l\H(r)\g(t))) 
drdp 5((l\P\g q (t))) 



(2nhr S((l\n Aa (xMt))) 



H(r)-((l\H(r)\ e (t)))) + 

n A M - ((M^AMleit))) 



and has all the properties of a projection operator: 
p(t)H(r) = H(r), p(t)p(t') 



P(t), 



p(t)n Aa ( r ,P) 



77 



Aa 



y,p), 



p(t))p(t) = 0. 



The obtained transport equations have the general meaning and can describe 
both weakly and strongly nonequilibrium processes of a quantum system with taking 
into consideration coupled states. In a low density quantum field Bose- or Fermi- 
system the influence of the average value of interaction energy is substantially smaller 
than the average kinetic energy, and coupled states between the particles are absent. 
In such a case the set of transport equations (|8.12|) , ( |S.13| ) is simplified. It transforms 
into a kinetic equation [ |27| in thermo field representation for the average value of 
the Klimontovich operator ((l\n(x)\g(t))): 



^«i|nfc(p)le(*)» = «i|A fc (p)|a 1 (t))> + 



T(t,t') 



Jdr Qltf)J n (grf)Q^tf))) b- g (p';t 



Using the projection operators method, this equation was obtained in [ |24 



In the next step we will construct such annihilation and creation superopera- 
tors, for which the quasiequilibrium thermo vacuum state vector is a vacuum state. 
Analysing the structure of quasiequilibrium statistical superoperator (p.8|), one can 
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mark out some part which would correspond to the system of noninteracting quan- 
tum A-particles. Let us write g q (t) in an evident form and separate terms which are 
connected with the interaction energy between the particles: 



ft,(t)=expj-$*(i)-y dr/3(r;t) 



x 



(8.20) 



E 



drdp f p 2 



2m A 



- / dr[3(r;t)H int (r) 



Using operator equality (A and B are some operators) 



i 

1 + JdT e r(A + B) Be -r(A+B) 



the relation for g q (t) can be rewritten in the following form: 



dr(3(r;t) / dr &(t)H int (r)(g q (t)y 



Q° q (t), 



$.21) 



where 



E 



X 



(8.22) 



dr dp 



P 
2m A 



n Aa ( x ) ~ ^Aa( x ;t)n Aa (x) 



or 



Q° q (t) = exp I ^(t) - J dr ftr; t) £ J j^b^x; t)n Aa {x) \ , (8.23) 



where b Aa (x;t) = 



2m A 



. Quasiequilibrium statistical su- 



peroperator g q (t) is bilinear on annihilation and creation superoperators a Aa (P) and 
a Aa (P), as well as on the non-perturbed part of Hamiltonian H . One can write the 
total quasiequilibrium superoperator as some non-perturbed part of g q (t) and the 
part which describes interaction of quantum particles in the quasiequilibrium state. 
Further, we introduce the following designation: 



Q q (t) = g° q (t) + g' q (t), 



where 



m = ~ j dr fart J dT gl(t)H mt (r)(g q (t)r T g° q (t). 



(8.24) 



(8.25) 
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Quasiequilibrium thermo vacuum states \g q (t))) and |£q(£))) are not vacuum states 
for annihilation and creation superoperators d Aa (P), d Aa (P), d Aa (P), d Aa (P). But 
for \g°(t))) one can construct new superoperators l Aa {P), ^f Aa (P), Ia^P), l Aa (P) 
as a linear combination of superoperators a Aa (P), d Aa (P) and a Aa (P), d Aa (P) in 
order to satisfy the conditions: 

7 AQ (P;t)K(t)» = o, «i|%U^;*) = o, 

W^;*)le!(*)» = o, ((MiL(P;t) = o. 1 ' } 

To achieve this let us consider an action of annihilation superoperators a Aa (P;t), 
d Aa (P;t) on quasiequilibrium state |^q(to))) : 

a Aa (P;t)\g° q (t ))) = f Aa (P;t - t )a+ a (P;t)\gl(t ))), 

d Aa (P;t)\g° q (t ))) = af Aa (P;t-t )d + Aa (P;t)\g q (t ))), (8.27) 

where superoperators d Aa (p;t), d\ a (p;t), d Aa (p;t), d Aa (p;t) are in the Heisenberg 
representation 

a Aa (P;t) = e-^d Aa (P)e^\ d Aa (P;t) = e"^ ~a Aa (P) e*** 

and satisfy commutation relations: 

[a Aa (P; t),d%p{P'\ t)] a = 5a,b5«MP ~ P') 
[d Aa (P; t), ~a + Bp {P'- 1)] a = 5 A>B S a>0 S(P - P') 
[*A a (P;t),~a B p(P';t)} a = [at a (P;t),a+ /3 (P';t)] a = 0. 

It is necessary to note that superoperators H(r), n Aa (x) are built on superoperators 
a Aa (p + §), a\ a (p - §), d Aa (p + f ), a\ a (p - §). Therefore, for convenience here a 
unit denotion was introduced for arguments like P = p ± |. This should be taken 
into account in further calculations where obvious expressions are needed. 

According to general relations of section 6 (|6.7p - (|6.19| ), we can introduce new 
operators l Aa {P)t), l Aa {P]t), n/ Aa (P;t), l Aa {P]t) via superoperators a Aa (P;t), 
a+ a (P;t),d Aa (P;t),d+ a (P;t): 

a (p t) - n ^*( p ;Mo) -+ x 



l Aa (P;t) = \ l + an Aa (P;t,t ) 



lta(P;t) = ^l + crn Aa (P;t,t )[d\ a (P;t)-ad Aa (P;t)]. (8.28) 
Relations (|8.28|) satisfy conditions (8.26). Here 

n Aa (p,q;t,t ) = n Aa (P;t,t ),= ({l\d+ a (P;t)d Aa (P;t)\gl(t ))) = 

((l|«Aa(P-f;*)«Aa(P+|;*)l^0)», 

is a quasiequilibrium distribution function of A-particle coupled states in momentum 
space p, q, which is calculated with the help of quasiequilibrium thermo vacuum 
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state vector f>q(tn))) Q.23|) . Function f Aa (P;t — in) in formulae (|8.27| ) is connected 
i Aa (P;t,t ) b 

f Aa (P;t-t ) 



with n Aa (P; t, in) by the relation 

n Aa (P;t,t ) 



1 + an Aa (P;t,t )' 

Superoperators 7 Aa (P;i) and 7 Ace (P;i), 7 AQ (P;t) and 7 Ace (P;i) satisfy the "cano- 
nical" commutation relations: 

liA a {P;t),iUP';t)] a = W q ^(p- p'), 

llA a {Pit),lw{P'it)] a = ?>A,B5 a ,p5{P - P'), (8.29) 

[iA a (P;t),i BP (P';t)] a = [iUP;t),iU p '^], = °- 

Inversed transformations to superoperators a Aa (P;i), a Aa (P;i) are easily obtained 
from (|8~28|): 



a Aa ( P ^) = \/ l + orn Aa(P^ t ^o, 



5+ a (P; t) = ^/l + an Aa (P;t,t ) [j+ a (P; i) + aj Aa (P; i)] . (8.30) 

7,4a (-^ i 0> 1ao.(P'i 0> 7y4a(-fj xL^i c °uld be defined as some operators of an- 
nihilation and creation of A-quasiparticle coupled states, for which quasiequilibrium 



thermo vacuum state |£q(io))) (|8.23|) is a vacuum state. In such a way, we obtained 



relations of dynamical reflection of superoperators a Aa (P;t), a Aa (P;t), a Act (P;i), 
a Aa (P;i) to new superoperators of "quasiparticles" J Aa {P',t), 7 Aa ,(P;i), 7 Aa (P;i), 

lAa(P-t). 

A set of transport equations ( fB.12j ), (|8.13| ) together with dynamical reflections 
( E>.28|) , (|8.30|) of superoperators in the thermo field space constitute the basis for a 



consistent description of the kinetics and hydrodynamics of a dense quantum system 
with strongly coupled states. Both strongly and weakly nonequilibrium processes of 
a nuclear matter can be investigated using this approach, in which the particle 
interaction is characterized by strongly coupled states, taking into account theirs 
nuclear nature [ 0, ^ [| . 

It is much sequential to describe investigations of kinetic and hydrodynamic 
processes of a nuclear matter on the basis of quark-gluon interaction. The quantum 
relativistic theory of kinetic and hydrodynamic processes has its own problems and 
experiences its impetuous formation [ [l], ^ |3|, ||], || . In the next section we consider 
one of the possible ways of describing the kinetics and hydrodynamics of a quark- 
gluon plasma. 

9. Thermo field transport equations for a quark-gluon plasma 

Investigation of the nonequilibrium properties of a quark-gluon plasma - QGP 
which can be created after ultrarelativistic collisions of heavy nuclei [ [j], §, [|. [|. 
|) 0' & S or l aser thermonuclear synthesis is topical from the point of view 
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based on the statistical approach ||, ||, |T], H ||, ||, H . Thus, 

it is important to construct kinetic and hydrodynamic equations for QGP because 
such a state of a nuclear matter is characterized by high temperature, large density 
and strong interactions between quarks and gluons which are described by chromo- 
dynamics [ [51], [52], |53| . At present there are certain achievements in this direction. 
The classical theory of transport processes in QGP, based on the relativistic Vlasov- 
Boltzmann equation, was proposed by Heinz [ [5||. On the basis of these works, 
transport coefficients for weakly nonequilibrium QGP were studied in [ 56, 57, 58 



was obtained in [ 59 



The Lenard-Balescu-type collision integral for the classical model of a quark plasma 

using the Klimontovich method [ |60], |6l[| . Hydrodynamics of 

The Vlasov-Boltzmann equations, like 



QGP is considered in papers [ f49l BO, [62], |63 



the quantum kinetic ones for Wigner distribution functions for quarks and gluons, 
were obtained in [ [19, Some interesting results were obtained in paper [ 



67J where the temperature behaviour of the kinetic coefficients of a gluon gas had 
been studied using the Green-Kubo formalism and <P 4 -model. 

However, it is necessary to note that the Vlasov-Boltzmann kinetic equation is 
correct for a rarefied plasma only (small densities). Thus, it is only the first step in 
the investigation of transport processes of QGP. 

In a dense high temperature quark-gluon plasma, which is characterized by strong 
interactions, kinetic and hydrodynamic processes are mutually connected and should 
be considered consistently. In this section we carry out a consistent description of 
the kinetics and hydrodynamics for QGP on the basis of a nonequilibrium thermo 
field dynamics using the method of nonequilibrium statistical operator 



, 34 



We will obtain generalized relativistic quantum transport equations of the consistent 
description of the kinetics and hydrodynamics for QGP. It should be also noted that 
problems of the description of nonequilibrium properties of QGP were considered in 
papers [ ||, |9|]. 

M 



Consider QGP with the Lagrangian from quantum hydrodynamics 
L = \f; v F; v + ^ + (i^ + | X a A a ^) 



(9.1; 



[IV 



duA a u + gf 



abc Ab ac 



where the fields of the matter are spinor quark fields & for which one employs the 
collective designation \P(x) with components where i = 1,2,3 (index of colour: 
red, green, yellow); / = 1, . . . , 6 (flavour indices: b, c, d, s, t, u); a = 1, . . . , 4 (spinor 
indices); are gauge vector fields (Jang-Mill's fields) that correspond to gluons; 
a = 1, . . . , 8 indices of colour; A a means eight Gell-Mann matrices satisfying the 
commutation relations 



A 61 
2 



if 



abc 



f 



abc 



are the structural constants of groups SU(3); g is the gauge constant of connec- 

o 



ct, x); 7^ are Dirac matrices [ |5T|, [52], ^3 



tion; dp = d/dx^, = (x° 

The nonequilibrium state of such a system is described by a relativistic quan- 
tum nonequilibrium statistical operator g(t) which satisfies the relativistic Liouville 
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equation admitting an obviously covariant form. Such a Liouville equation for a 
thermal quantum field system in a covariant form in the interaction representation 
was written down in [ 57, 7C\ on the basis of the Tomonaga-Schwinger equation [ 



7T|, |72| in the next form: 



-Qmt(<r) -i[#int(a:),ft„t((x)] = 0, (9.2) 



8a(x) ih 

where the nonequilibrium statistical operator g(a) is defined on an arbitrary space- 
like surface a(x) [ [71], [7^, [73|]. When the surface a(x) tends to the plane t = const, 
equation ( |9.2D transforms into the quantum Liouville equation in the interaction 
representation, 

d 1 

Qj.&rt(t) ~ ^ [H int (x), g- mt {t)\ = 0, 

where g- m t{t) is a nonequilibrium statistical operator in the interaction representa- 
tion, given on the plane t = const: 

H int (t) = exp | _ r^^| Ant exp jr^^j , Qmt(t) = exp | _ r^^| g(t), 

Hq and H mt are noninteracting and interacting parts, correspondingly, of the total 
Hamiltonian H of the system. The Hamiltonian of a quark-gluon system correspond- 
ing to the Lagrangian ( |9.1| ) was obtained in [ [74[] using the Coulomb gauge. We will 
represent it in the following form: 

H = H + H int , (9.3) 
H = l -P?P? + \d n A^d n Al 

flint = gr hc diA a n A\A c n + ^f^f^AlA^AfAl - ]-A a AA a + 



(iflk + |AMf) 7 z + m - |AMS7o 



selecting the "free field" Hamiltonian Hq and the one which describes an interaction 
between quark and gluon fields. Here m is a colour independent mass matrix for 
quark flavour indices, means a canonical momentum conjugated to gluon field 
Af, while zeroth components Aq, like in electrodynamics, are not independent and 
should be inserted in H int after the solution of the equation of motion: 

AA a = 9 r hc A\ (P t c + d t A c ) + !<^ + 7oA^. 

On the contrary, proceeding from electrodynamics, it is possible only in the form 
of infinite series, so that Hamiltonian (|9.3|) really consists of an infinite number of 
vertices. 

To solve the relativistic Liouville equation ( |9.2|) , boundary conditions should 
be set. Let us search the solutions which depend on time through a certain set 
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of observable quantities only, the number of which is sufficient for the descrip- 
tion of the nonequilibrium state of a system, using the method of a nonequilib- 
rium statistical operator. For this purpose, we introduce the infinitesimal source 

— v(Qint( a ') — f?q nt (°")) i n the right-hand side of equation ( |9.2|) , which corresponds 

to the boundary condition ft nt (cr) — > f? q nt (cr) with a — > -co and, according to the 
formalism of nonequilibrium thermo field dynamics, write down this equation in the 
thermo field representation: 

' IM<t)» - ^H- mt (x)\ eint (a))) = - V (\Qint(cr))) - |^ nt (a))) ), (9.4) 



Sa(x) ih 

where the source selects retarded solutions with respect to a reduced description of 
the nonequilibrium state of a system. |^ q nt (cr))) is a quasi-equilibrium state- vector; 
Hi n t(x) is the Hamiltonian of interaction between the quark and gluon superfields, 



H Q (x) = H (x) - H Q (x), H- mt (x) = H- mt (x) - H n 



where Hq and H bilinearly depend on Bose superfields and A ^ (without and with 
a tilde line) with corresponding commutation relations analogous to those for the 
Bose fields [ ^TJ; H int and H int are Hamiltonians composed of gluon Bose superfields 
A a ^ and A a , of quark Fermi superfields \P + and \I/ without a tilde line, and \P + and \P 
with a tilde with corresponding commutation alignments which are similar for both 
Bose A a ^ and Fermi fields W [ 51 1. We shall write "Schrodinger" equation (|9.4j ) in an 



integral form, introducing 

J dV exp( V Q a/(7 )T(a,a') (-^ - ^int(*')) \g^\a'))), 

— oo 

where Q a ',a means the volume enclosed between the surface of integration a'(x') and 
the second one a(x); T(a,a') is an evolution operator, 




T(cr, a') = exp + 

The quasiequilibrium statistical superoperator g q (cr) in \Q q (cr))) = g q (a)\l)) is de- 



fined generally [ |27|, |33|, |34| from the entropy extreme with keeping the normal- 
ization and under the conditions that observable quantities (p n )' = ((l\Pn\ (>(&))) 
are given. There are the following characteristic values for QGP: mean density of 
the quark barion charge {b fl (x)) t , mean densities of the colour currents (j^ 
for quarks and (j^ J X )Y f° r gluons and also values of the densities of the energy- 
momentum tensor (T^ u (x)Y for quarks and (T^(x))' for gluons, in which densities 
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+ satisfy the local conservation laws, 



0. 



d, K „(*) + # „(*)] + ^"^(x) K „(*) + J g c i M (x)] = 



0„i£,(s) + 5 r 6c 4(x)F^(x) - 5 (j« „(*) + 4 „(x 



x) 



(9.6) 



0. 



The densities of operators b^(x), j* (x) + A x )i T^ v (x) + T&(x) are "slowly- 
changing" and their mean values on the long time scale satisfy the relativistic equa- 
tions of hydrodynamics; it is necessary to co-ordinate with them the kinetic equations 
for quark and gluon distribution functions. Local laws ( |9.6|) impose some restrictions 
on the kinetic processes, and what is more, their role is considerably important at 
large densities and strong interactions. It shows that the kinetics and hydrodynamics 
of a quark-gluon system are strongly correlated. That is why, in writing down the 
kinetic equations for such systems, it is natural to choose the reduced description 
of the nonequilibrium state in such a way that the proper dynamics of conserved 
quantities is taken into account automatically. 

We introduce Wigner operators for quarks and gluons [ |6B], [7S| to obtain 

the kinetic equations: 



f(x;p) 



T— 774 eX P { —PV }U(X,X_)& (x_) $ + (X + ) U (x+, X) 

(27m) 



1 

m 1 



is the Wigner operator of quark density; 



G^(x;p) 



d 4 y 



1 



U(x,x-)F lxX (x-)U(x-,x 



(9.7) 



{9.1 



U(x, x + )F^(x + )U(x + , x) 



is the Wigner operator which is connected with the density operator of an energy- 
momentum tensor for gluons: 



-gV Sp G$(x)). 



d p (SpG^y 

Here U(x', x) is an operator which was introduced in paper 



(9.9) 




dz» AJz 




(9.10) 



A;(x 



U(x', x) is a connection operator, where the integration path z(s) = x + s(x' + x), ^ 
s ^ 1 lies in the plane t = const. The densities of barion charge bJx), quark colour 
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current (x) and energy- momentum tensor Tf} u (x) are defined via the Wigner 
operator of quark density, for example, 



(9.11) 
(9.12) 
(9.13) 



b^x) = dp (Sp7 M /(a;;p)), 



X 



h 

d A p (Sp-A 7 M /(x;p)), 
d 4 p (Sp7 M /(a;;p))P I/ . 



To provide a consistent description of a quark-gluon system, we define the quasi- 
equilibrium statistical superoperator g q (cr), g q (a)\l)) = |f? q (cr))) in a standard way 
with keeping to the normalization condition and the following requirements of the 
averages [ 571, R3, M: 



((lK(x)\g(cr))), 
({l\% v (x)\g(a))), 

Then g q (cr) takes the form: 
g q (a) = exp {-S(cr 

where 



({l\f(x;p)\g(°))), 
((l\G, v (x;p)\g(a))). 



S(a) = + J da» [p%„ - + 



(9.14) 



(9.15) 



(9.16) 



da 11 I dp a(x;p)f(x;p)+uj u (x;p)G flu (x;p) 



cr(x) is an arbitrary space-like surface passing through point x: dcr^ = dan^ is a sur- 
face element vector; is a normal vector (da = d 3 x); P v (x), £,(x), a(x;p), cu u (x;p) 
are the Lagrange factors that are defined from the self-consistency conditions: 



({l\K(xMa))) 
((l\f^x)\g(a))) 
((l\JXx;p)\g(a))) 
((l\G^(x;p)\g(a))) 



{{l\b,(x)\g q (a))), 
((l\T, u (x)\g q (a))), 
((l\f(x;p)\g q (a)}), 
((l\G^(x;p)\g q (a)}}, 



(9.17) 



= (3u IJL1 (3 = 1/T(x) denotes inverse invariant temperature, u^(x) is a local hy- 
drodynamic velocity; £ = (3fi, /i a (x) is a local chemical potential of a-sort quarks 
(a: b, c, d, s, t, u). Parameters a(x;p) and u u (x;p) are conjugated to averages 
((l\f(x;p)\g(a))), ((l|G AtJ/ (x;p)|^(a))). T^Qc ), 6 M (x), f(x;p), G^(x;p) are superop- 
erators constructed according to d9.7|) -( p^l~3|) and the thermo field dynamics formal- 
ism [[27|, on Bose superfields A a (x) and Fermi superfields <? r+ , \P. <P + (cx) is calculated 
from the normalization condition: ((l\g^(a))) equals to 1. Accordingly, superopera- 
tor g+ (cr) is built on operators T M „(x), b^(x), f(x;p), G^ v (x;p) which are expressed 
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via Bose superoperators A'tix) and Fermi superoperators \P (with a tilda line). 
Taking into consideration the structure of quasiequilibrium statistical superopera- 
tor (|9.15 ) and ( |9.16 ), one can write the nonequilibrium thermo vacuum state vector 
31) in the interaction representation in the following form: 



lfet(a))) = |^ nt (a))) - / dVe^r(a,a')x 



5 



1 ^. 



(9.18) 



5a'(x') ih 

Using the rule [§§ ^ f a da»(x') A„(x') = ^A^x), for (5S in V))/(<V« 

d r 



we obtain 
5 



5a' (x' 



CJ"int I _l\ 



(J 



dx', 



P»T*(x')-tb™\x') + 



d 4 p [a(x';p)r t (x';p)+^(x';p)G^(x';p)] 



The action of the operator ——Hi nt (x) on S mt (a') reads: 



-H int (x)S mt (a') 
in 



da^(x") P u f^(x") - &™\x") 



+ 



da^x") I d 4 p a(x"; p)r\x" ] p) + oo»(x"; p)G^(x"; p) 



where 



1 



1 



ftfV) = [H mt (x'),V»\x")} , r\x" ]P ) = -- [H int (x'),r\x";p)} , 
T™Hx") = -1 [H mt (x'),T^(x")} ,G%(x";p) = -1 [F int (*'), G^(x"; p)] . 



ih 



Using the nonequilibrium thermo vacuum state vector (|9.18|) , we obtain a system of 
coupled equations for averages ( 9. 14j ) in the interaction representation: 







(9.19) 



P ^((lK(x)\g- mt (a))) = 0, 

P^((MW\QUa)}} = P » 1 ir(( 1 \WK t (°))} + (9-20) 



dx. 



d 



p»— I d 4 x'e^((l\f, u (x)T(a,a')\x 



i 

/ d ^ e - rSmt(CT,) {^y - i H ^} s m VK 5int(CT VV)>>, 
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p ■ D(x)((l\f(x }P )\g int (a)}} = p ■ D(x)((l\f(x,p)\g^(a)}} + (9.21) 
p-D{x) [ dVe' ?n -'-((l|/(x,p)T(or,a')| x 



1 

o 

p ■ D(x)((1\G^(x, P )\qU<7))) = P ■ D(x)mG^x )V )\Qf(a))) + (9.22) 

(7 

p-D(x) J dV e^((l|G^(x,p)T(a,a , )| x 



o 

while the space-like surface cr(x) tends to the plane t = const. Here D(x) and D(x) 
are covariant derivatives. They act on the Wigner operators of quarks and gluons, 
correspondingly: 

D lx {x)f{x,p) = + hhcA^x^j f(x,p), 

where A^{x) = A a ^(x)T a with (T a ) bc = —ihf abc is an 8 x 8 matrix. The tildian 
covariant derivative acts on the Wigner operator of gluons 

d 

D a (x)G fll/ (x,p) = Q^r G ^( x 'P)jfc l A »,G^ u (x,p)} 



with a commutator between two 8x8 matrices {A a ) mn = A a a (T a ) mn and (g^j . An 

action of these covariant derivatives on the Wigner operator of quarks p- D(x)f(x,p) 
and on the Wigner operator of gluons p ■ D(x)G ^(x,p) was calculated in papers [ 
6^1 , The set of equations (|9.19|) - (]5.22|) without taking into account the barion 



charge transport ((l\b fl (x)\g(t))) and the total momentum tensor ({l\T tM (x)\g(t))) 
turns to a connected system of kinetic equations for the average values of Wigner 
operators of quarks ((l\f(x,p)\g(t))), gluons ((l\G^(x,p)\g(t))) . It generalizes the 
results of papers [ |63], |65], |66| . 

We have obtained a system of coupled relativistic transport equations of a con- 
sistent description of the kinetics and hydrodynamics for QGP in thermo field rep- 
resentation. This system of equations is strongly nonlinear and could be used to 
describe both strongly and weakly nonequilibrium states of a system. However, it 
should be noted that the following transformations for the use of the system of 
relativistic transport equations (|9.19| )- (|9.22| ) have to be made. Since the quasiequi- 



librium thermo vacuum state vector |f?q(<r))), which is used in averaging equations 
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fl9.19|) -( |9.22| ), is not a vacuum state vector for Bose superfields A®, A® of gluons 

and Fermi superfields & + , \P, \P + of quarks, it is necessary to construct Bose 
superfields of gluons and Fermi superfields of quarks for which it is a vacuum state 
vector, as in paper [ |27|. Such superfields will depend on the parameters of state 
P u (x), (3(1, a(x;p), u u (x;p), besides, all the superoperators T^(a;), b^(x), f(x;p), 
G^(x;p), Hi nt (x), T(a,a') in a system of relativistic transport equations (|9.19|) - 
( |9.22D should be expressed in a such way where \g q (cr))) is the vacuum state vector. 
Besides, the use of space-like surface a(x) is justified at the calculation of generalized 
equations in the invariant form. Since the final results do not depend on the choice 
of surface a, one has to rewrite equations ( |9.19| )-( |9.22| ) on the surface directed to 
the plane t = const. So, it is useful to choose r\ = eV' 1 , where V denotes the volume 
occupied by system [ [7(J. Then VL a i a = (t' — t)V and r/fV'o- = s(t' — t). Thus, in 

transport equations (|9.19|) - (|9.22|) , J d A x' exp ( nVL^^ ) — > J dt' exp [e(f - t)) 

— CO — oo 

and superoperators ft n t(t), Q^it), T(t,t') will be given on the plane t = const. Such 
questions will be considered in the next paper, in particular, in the investigation of 
the relativistic transport equations for weakly nonequilibrium QGP. 



10. Conclusions 

In this paper the conception of nonequilibrium thermo field dynamics on the ba- 
sis of a nonequilibrium statistical operator has been applied to the construction of 
transport equations of dense quantum systems. Hydrodynamic equations in thermo 
field representation have been obtained for both strongly and weakly nonequilib- 
rium processes. Transport cores, connected with transport coefficients of viscosity 
and thermal conductivity, have been defined, too. They are calculated with the 
help of a quasiequilibrium thermo vacuum state vector. The last one is a vacuum 
state vector for annihilation and creation superoperators which depend on thermo- 
dynamic parameters F n (t) of a system. The consistency of both annihilation and 
creation superoperators and the corresponding vacuum state is accomplished then 
by the mentioned approach. In particular, it is achieved by the generalization of 
the Bogolubov transformations in thermo field dynamics [ [32|] for a nonequilibrium 
case. General transport equations within nonequilibrium thermo field dynamics - 
nonequilibrium statistical operator made it possible to obtain an equation of a gen- 
eralized description of the kinetics and hydrodynamics for a dense quantum system 
with strongly coupled states. In such a case, fAa(x;t) - the Wigner function of 
the A-particle coupled state, and the average value of the total energy density op- 
erator {{l\H(r)\g(t))) have been chosen as parameters of a shortened description. 
Particle annihilation and creation superoperators for the quasiequilibrium thermo 
vacuum state vector of noninteracting particles were defined in (|3.22|) . It gives us 
the possibility for the construction of a diagram technique in the calculation of the 
corresponding transport cores. These investigations and calculations are important 
for actual nuclear systems 
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We have considered a method for obtaining generalized transport equations for 
QGP - one of the nuclear matter states. These equations were obtained in the most 
general form. To be used, their structure needs a lot of transformations, especially 
for transport cores. The consistent description of the kinetics and hydrodynamics 
is based on a set of fundamental parameters of a shortened description: the average 
values of Wigner operators for quarks and gluons ((l\f(x)\g(t))), {{1\G ^ u (x)\g(t))) 
and the average values of the density operator of barion charge ((l|6 M (x)|^(t))) and 
the total energy-momentum tensor ((l|T Atl/ (o;)|^(t))) of quarks and gluons which 
constitute the basis for the hydrodynamic description. The set of equations obtained 
permits the investigations of weakly nonequilibrium processes and kinetic equations 
like the Boltzmann-Vlasov or Lenard-Balescu ones for diluted QGP. 

The problem of the investigation of transport coefficients: viscosity and thermal 
conductivity as well as excitations of QGP still remains. It might be considered in 
view of the obtained transport equations. This will be the subject of our future 
work. 
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